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ABSTRACT 

The Bar is the most productive region of the Small Magellanic Cloud in terms of star formation 
but also the least studied one. In this paper we investigate the star formation history of two 
fields located in the SW and in the NE portion of the Bar using two independent and well tested 
procedures applied to the color-magnitude diagrams of their stellar populations resolved by means 
of deep HST photometry. We find that the Bar experienced a negligible star formation activity in 
the first few Gyr, followed by a dramatic enhancement from 6 to 4 Gyr ago and a nearly constant 
activity since then. The two examined fields differ both in the rate of star formation and in the 
ratio of recent over past activity, but share the very low level of initial activity and its sudden 
increase around 5 Gyr ago. The striking similarity between the timing of the enhancement and 
the timing of the major episode in the Large Magellanic Cloud is suggestive of a close encounter 
triggering star formation. 
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1. INTRODUCTION 

The Small Magellanic Cloud (SMC) is a funda- 
mental laboratory to study the evolution of dwarf 
irregular galaxies (dirr's). The SMC is the clos- 
est member of this class of systems, has a current 
metallicity (Z ~ 0.004 as derived from HII regions 
and young stars) similar to that of the majority of 
dirr's and a mass (between 1 and 5 x 10^ Mq, e.g. 
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Kallivavalil et al. 20061 and references therein) at 
the upper limit of their range. Moreover the SMC 
is a member of a triple system, a circumstance 
that favors studying the modulation of the star 
formation activities driven by interactions. 

In order to derive the detailed, spatially- 
resolved star formation history (SFH) of the SMC 
we have started an international long-term project 
to study t he evolution of t he SMC in space and 
time (see Tosi et al. 2008f) . Our strategy is to 



achieve high photometric depth and spectroscopic 
resolving power over a large field of view by com- 
bining datasets from the ground and space. We 
are using the Hubble Space Telescope (HST), the 
Very Large Telescope (VLT), and the VLT Survey 
Telescope (VST) to observe a large sample of field 
stars and clusters across the SMC. These data will 
allow us to constrain the global SFH as well as the 
existence of chemical and age gradients. 

For the cluster analysis, we have already pre- 
sented deep photometry with HST's Advanced 
Camera for Surveys (ACS) of seven intermediate- 
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age and old populous clusters (jGlatt et al ■l2008allbl . 
20091 l201l[ ). In combination with our VLT data 



Chiosi fc Vallenaril (|2007^ derived the SFH in the 
vicinit y of a few SMC cl u sters. Finally, et al 



we find a complex age-metallicity relation for these 
clusters with a considera ble spread in metal licity 
at any given age (see e.g. Glatt et al]|2008bl ). 

Concerning the SMC field analysis, our plan 
is to have Color-Magnitude diagrams (CMDs) 
several magnitudes fainter than the oldest main- 
sequence (MS) turn-off (TO) for the entire galaxy. 
To this pu rpose we have obs erved six fields with 
the ACS (ISabbi et al.l l2009l) . sampling regions 
characterized by different stellar and gas densi- 
ties in the SMC Bar, in the Wing in the direction 
of the LMC, and in the outskirts (see Fig. [1]). A 
preliminary SFH analysis of suc h fields has been 



presented in ISabbi et al. (I2OO9I) . We also have 



an ongoing Guarant eed Time O bservation pro- 
gram at the VST ( Ripepi et al.l [2OO6 ) designed 
to cover with deep photometry the whole SMC 
and the Bridge connecting it to the LMC. These 
CMDs will allow us to infer for the first time the 
SFH of the whole SMC over the entire Hubble 
time, covering a much larger ar ea with con sid- 
erably better image quality than IZaritskv et al.f s 
(120021 ) data. We will derive the SFHs from the 
CM Ds using the synthet ic CMD technique (see 
e.g. lTolstov et al.lboOQ: Cignoni fc To si 2010,, and 
references therein). SFHs of some SMC fields 
have been derived and presented by other au- 
thors, based on ground-based observations or HST 
studies of small individu al regions (see Fig. [T|). 



Harris fc Zaritskvl ()2004l) derived the SFH of the 



SMC over 4° x 4.5° to a depth of V < 21 using the 
Magellanic Cl oud Photomet r ic Su rvey (MCPS) 
UBVI catalog IZaritskv et al.l (|2002[ ). This is cur- 
rently the most spatially extended study of the 
galaxy, but does not reach the oldest MSTO. The 
most comprehensive study of the old populat ion of 



the SM C to date was carried out by Haschke et al 



(|2012[ ) using RR Lyrae stars from the Optical 
Gravitational Lens ing Experiment (OGLE-HI; 
lUdalski et al.|[2008[) . They find a uniform metal- 
licity distribution across the SMC with a sp r ead o f 
more than 1 dex in [Fe/H]. lOolphin et all (|200ll ) 
analyzed the stellar content of the SMC halo, in 
a region close to the globular cluster NGC 121, 
using both HST Wide Field Planetary Camera 2 
(WFPC 2) and ground based d ata. Again with 
WFPC2 iMcCumber et al" Idlooi) studied the stel- 
lar content of a small portion of the SMC Wing. 



(|2007l ) and iNoel et al.l (|2009l ) presented a deep 
ground-based study of 12 fields of the SMC, avoid- 
ing the densest regions, because of their high 
crowding conditions. 

In this paper we present the SFH of SFHl and 
SFH4, the two most central fields of the six SMC 



regio ns we observed with HST/ ACS (Sabbi et al 



20091) . The apparent distances from the SMC op- 
tical center are about 24' and 1° 52' respectively. 
SFHl is located in the SW portion of the SMC 
Bar, where the stellar density, gas and dust con- 
tents are highest, while SFH4 is located to the NE 
of the SMC center, at 24' south of NGC 346, the 
most active star- forming region in the SMC. 

For a better assessment of the intrinsic theo- 
retical uncertainties, the SFH is derived using two 
completely independent procedures for the appli- 
cation of the synthetic CMD method. We com- 
pare the two methods here and discuss the corre- 
sponding results. The other ACS fields o bserved 
by us and described by ISabb i et al.l (120091 ) are be- 
ing treated in the same way, and their SFH will 
be presented in a forthcoming paper. 

We briefly describe our data in Section 2. The 
two procedures for the SFH derivation are sum- 
marized in Section 3, together with the results of 
their application to SFHl and SFH4. Similarities 
and differences between the resulting SFHs are dis- 
cussed in Section 4, while in Section 5 we compare 
our findings with previous literature. Concluding 
remarks follow in Section 6. 

2. Photometry and CMDs of SFHl and 
SFH4 

The data for our six SMC fields were ac- 
quired with the ACS Wide Field Channel between 
November 2005 and January 2006 (GO-10396; P.I. 
GaUagher) with the F555W and F814W fihers. 
The data reduction was perfor med with the pro- 
gram img 2xym_WFC.09xlO (|Anderson fc Kind 
20061 ). and the resulting magnitudes were cali- 
brated in the Vegam ag photometric system using 
Sirianni et al.l (|2005l ) recipes. For sake of simplic- 
ity, from now on we will refer to the tofsssw E^nd 
TTipsiiW magnitudes calibrated in the Vegamag 
system as V and I, respectively. 

Extensive artificial star experiments were per- 
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Fig. 1. — Spatial dis tribution of the six observed fields (red s ymbols) together w ith the observatio ns from 
iDolphin et aD (l200ll) Cblue symbols ), iMcCumber et all (|2005[ ) (magenta symbol), iNoel et all (|2007t ) (green 
symbols). IChiosi fc Vallenaril ( 2007f ) (yellow symbols), superimposed on the DSS image of the SMC. 
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formed to test the level of completeness and the 
photometric errors of the data. They followe d the 
approach described in I Anderson et al. ( 20081 ). and 
the artificial stars were searched for with exactly 
the same procedure adopted for the real stars. We 
considered an artificial star as recovered if its in- 
put and output positions agree to within 0.5 pixels 
and the fluxes agree to within 0.75 mag. As done 
in the photometric analysis we also required that 
each star was found in at least three exposures 
with a positional error < 0.1 pixels per filter. 

Details on both the data re duction and the art i- 
ficial star tests can be found in ISabbi et al.l(|2009f ). 

The final SFHl and SFH4 catalogs contain 
about 29200 and 17300 stars respectively, and the 
corresponding CMDs are shown in Figl2l where 
the photometric errors are also plotted These 
CMDs reach almost four magnitudes fainter than 
the oldest MSTO, thus allowing to study the evo- 
lution of the regions over the whole Hubble time. 
The sub-giant branch (SGB) is well populated and 
its brightness extension, much larger than the pho- 
tometric error at that magnitude level, suggests a 
prolonged star formation between 3 and 12 Gyr 
ago. The red giant branch (RGB) and the red 
clump (RC) also are well populated. In addition to 
these intermediate-age and old components, both 
CMDs show a MS blue plume and a blue loop (BL) 
sequence typical of late-type dwarf galaxies, cor- 
responding to young high- and intermediate-mass 
stars. 

Finally, the CMDs do not show any signifi- 
cant population of stars on the right of the lower 
MS, the CMD locus occupied by pre-main se- 
quence (PMS) stars. This, in turn, suggests a 
negligible activity in the last 50 Myr (the aver- 
age time spent by a sq lar-like star in PMS; see, 
e.g.. lCignoni et al.ll2010l ) compared to more active 
regions of the SMC l i ke NGC 346 or NGC 602 

20091: 



see e.g. Nqta et al.l 2006 : Cignoni et al 



Carlson et al.ll201lUCignoni et al.ll2011l) 



^To be conservative, the plotted error bars correspond at 
each magnitude level to the larger value of the error result- 
ing from the photometric package and from the artificial 
star tests. 



3. SFH derivation: the two synthetic 
CMD approaches 

The SFH of SFHl and SFH4 has been de- 
rived with the synth etic CMD m e thod. As 
recently reviewed by Tolstoy et al. ( 2009t ) and 
Cignoni fc Tosi I (|2010[ ). this approach was first 
appli ed to nearby res olved galaxies twenty years 
ago (jTosi et al.l Il99l[ ) and is now recognized as 
a powerful tool to disentangle the complex his- 
tory of resolved galaxies. To both fields we have 
applied it f ollowing two inde pendent procedures: 
Cole's (e.g. ICole et al. 2007) and a combination 
of Cignoni's (e.g. ICignoni et al.M2006.) and the 
Bologna (e.g.lTqsi et al .lll99ll:lGre ggio et al.l ll998t 



Angeretti et al. 2005[) codes. These tecniques su 
persede the classical isochrone fitting approach, 
allowing us to explore a much wider parameter 
space and to incorporate statistical and obser- 
vational uncertainties. Nonetheless their use is 
limited by the reliability of the adopted stellar 
tracks and color transformations, as well as the 
nature of the IMF (e.g., §4.1) and quality of the 
data. The recovered history is the best of all pos- 
sibilities, but this does not necessarily imply that 
it is the actual solution nor that it is unique. 

The procedures used here have already been 
tested in comparisons with each other and similar 
methods. Both groups partic ipated in the Coim- 
bra experiment in 2001 (see ISkillman fc Gallart 
20021 and references therein), Cignoni's method 
has been applied to the derivation of the Solar 
Neigh borhood SFH (|Cignoni et al.ll2006t Icignoni 
20061 ) and in combination with the Bologna 
code to t he analysis of the SMC y oung clusters 
NGC 346 (ICignoni et al.l2010l l201ll) and NGC 602 



( Cignoni et al. 2009f) . Cole's method has been 



compared with two different cod es in the anal 



ysis o f the dwarf irregular IC 1613 (ISkillman et al 



2003 ) and the Cetus dwarf spheroidal (|Monelli et al, 
20101 ). Since the procedures are in continuous evo- 
lution, to improve both the reliability of the results 
and the efficiency of the computations, it is useful 
to perform further comparisons on deep and tight 
CMDs such as those of SFHl and SFH4. 

In the following we briefly summarize common- 
alities and differences of the two methods. In 
both cases the synthetic CMDs have been built 
using the results of the artificial star tests men- 
tioned above to assign photometric errors and in- 
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Fig. 2.— CMDs of the SFHl (left panel) and SFH4 (right panel) field observed with the ACS/WFC. The 
solid and the dashed red lines indicate the 50% and 25% levels of completeness, respectively. Formal errors 
on the estimated photometry are shown on the left side of each CMD (see text for details). 
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completeness to the synthetic stars. 

3.1. Bologna procedure 

The Bologna approach involves the compari- 
son of the observed CMD with a library of syn- 
thetic CMDs computed with different values of 
metallicity, initial mass function (IMF), binary 
fraction, distance modulus and reddening. Mod- 
els have been calculated with the latest Padova 
stella r models ( Bertelli et al. I l2008i iBertelh et al 



20091 ) for masses between the hydrogen burn 
ing limit and 20 Mq. Theoretical temperature 
and luminosity are transformed to the observa- 
tional plane by means of t he relations obtained by 
Origlia fc Leithereil (|2000f ) for the HST Vegamag 
photometric system. The synthetic CMD is cre- 
ated following a classical Monte Carlo procedure: 
1) stellar masses and ages are randomly extracted 
from a time independent IMF and a star formation 
law; 2) the stellar tracks are interpolated deriving 
the absolute photometry for the synthetic pop- 
ulation; 3) corrections for the distance modulus 
and the foreground reddening are applied. Then 
the synthetic CMD is degraded to match both 
the completeness profile and the photometric er- 
ror propert ies as derived from extensive artificial 
star tests ( Sabbi et al. 20091 ). To be conserva- 
tive, we have limited our fitting procedure to stars 
brighter than V — 23, whose completeness is bet- 
ter than 40% and photometric errors smaller than 
0.05 mag. 

The full SFH of a galaxy can be a complex func- 
tion of time. To make the problem manageable, we 
were forced to limit the range of parameter space 
that our models cover on the basis of previous re- 
sults and indications from our data. In order to 
reduce the computational time, the star formation 
rate is parametrized as a set of constant values 
over adjacent temporal steps: a generic CMD is a 
linear combination of chosen basis CMDs, where 
each basis is a Monte Carlo extraction from a step 
star formation episode. The duration of each step 
is chosen in relation to the evolutionary timescale 
of the average stellar mass of the step, ranging 
from 50 Myr (approximately the MS timescale for 
an 8 Mq star) to 1 Gyr (as representative of the 
theoretical precision at 10 Gyr, i.e. 10%) for ages 
above 3 Gyr. The star formation rate is constant 
within each step. In this work we deal with 90 
basis CMDs, each representing the synthetic pho- 



tometry of one of 18 age bins and one of five 
metallicities {Z = 0.008, Z = 0.004, Z = 0.002, 
Z = 0.001, Z = 0.0004). In order to reduce the 
Poisson fluctuations, all partial CMDs are gener- 
ated with more than 10^ solar masses. 

Within the framework of the adopted stellar 
tracks and atmosphere models, the most likely so- 
lution to the underlying SFH is the one which min- 
imizes the differences between data and synthetic 
star counts over strategic regions of the CMD. 
The degree of likelihood is assessed through a 
function of the residuals. Our experience suggests 
that the performance of such minimization is very 
sensitive to the adopted CMD binning scheme. 
Both fine and coarse grids o ffer advantages, as we ll 



Cignoni et adbOOel) . 



as disadvantages (see e.g 
Along the MS, a fine grid is mandatory to study 
the old stellar generations, since they are tightly 
packed together in the CMD. However, such a so- 
lution would pay the penalty of under-weighting 
the star counts along the upper MS which is Pois- 
son dominated. Vice versa, a coarse grid would 
be more adequate to map the recent activity, but 
it would allow a worse resolution at early times. 
The situation is even worse with the evolved stars, 
where the theoretical uncertainties are typically 
larger than for MS stars: the CMD position of 
horizontal branch (HB) and RC stars is affected 
by a complex interplay of age, metallicity , mass 
loss and helium content (see e.g. [Castellan i et al. 
l2000t ): not only the color but also the shape of 
the RGB are strongly affected by the color trans- 
formations from the theoretical to the observa- 
tional plane and by the mixing length parameter 
(in a way which depends on the used wavelength); 
the BL morphology is affected by the He burning 
cross sections (especially the ^^C(a,7)^^0) as well 
as by the efficiencies of external convection, over- 
shooting and mass loss. Nevertheless, the study of 
evolved stars can prevent the ambiguities in mod- 
els restricted to MS star counts, for example, when 
photometric errors and incompleteness hinder the 
possibility to trace the old activity by means of low 
mass MS stars. Given these issues, the problem of 
the grid choice can be solved only by using a vari- 
able grid spacing. Several tests were conducted for 
finding the optimal configuration. Figure |3] shows 
our best scheme. To capture all the information 
contained along the MS, the grid spacing shrinks 
with luminosity, balancing the longer evolutionary 
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Fig. 3. — CMD grid used to derive the SFH via the Bologna procedure. 



times of lower mass stars and, at the same time, 
providing sufficiently good statistics in the upper 
MS. Towards the red part of the CMD the average 
grid spacing is coarser: this helps to take the un- 
certainties into account, while preserving the num- 
ber of stars in specific phases. 

Finally, the is minimized by means of a 
downhill simplex algorithm. In order to avoid local 
minima, the simplex is re-started from thousands 
of initial random positions. A bootstrap method 
is used to assess the effect of random errors. The 
search of the best SFH is repeated for each boot- 
strapped data set, producing a distribution of best 
solutions. The errorbars on the final SFH repre- 
sent one standard deviation using 100 bootstraps. 

With this procedure we have analyzed the CMD 
of SFHl & SFH4. As a first step we have de- 
rived the SFH assumin g a two-exponent power law 
IMF rsee lKroupal200ll) with exponent s=2.3, close 
to Salpeter's 195^ for masses above 0.6 and 
s = 1.3 below, and a binary fraction (with primary 
and secondary mass extracted from the same IMF) 
of 30% (this parameter has no substantial impact 



on the fit quality). To reduce the parameter space, 
only three metahicities, Z = 0.004, Z = 0.002 and 
Z — 0.001, have been adopted. The distance mod- 
ulus and the foreground reddening were allowed to 
vary freely. The solutions with the lowest reduced 

are listed as SFHl-A and SFH4-A for the re- 
spective fields. 

In order to test how much the details of the 
solution depend on the allowed metallicity range, 
we calculated a second solution where the Padova 
tracks with Z = 0.0004 and Z = 0.008 were also in- 
cluded. These solutions are presented as SFHl-B 
and SFH4-B. These are generally consistent with 
the solutions obtained with the restricted metal- 
licity range but provide a somewhat better match 
to the observed CMD as described in ii4.1lbelow. 



3.2. 



Simulated annealing procedure (Cole's 
method) 



^Salpeter's IMF has exponent s=2.35 and the form 
Jq^°° m(p{m) dm = 1 and (f>(m) = in~^. 



The code developed by Cole has many features 
in common with the Bologna code, but also incor- 
porates some differences in the treatment of the 
Hess diagram data and in the calculation of the 
merit function that measures the relative likeli- 
hood of various solutions. Because the errors in 
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extracting detailed information from a CMD are 
dominated by systematic effects, the application 
of multiple fit procedures to the same data can 
provide important insight into which features of 
the SFH are robust and which may be artifacts. 

Cole's SFH-fitting code begins with a set of the- 
oretical isochrones interpolated to a fine grid of 
age and metallicity in order to create a synthetic 
CMD with no gaps. The most recent isochrones 
from the Padova models are used; all calculations 
are made in the HST Vegamag system to minimize 
the possibility of introducing errors in transforma- 
tion equations. The synthetic CMDs are binned 
in age and metallicity (Z) to increase the speed of 
computing and to avoid "overfitting" noise in the 
CMDs. We begin with bins that are evenly spaced 
by wO.lO in Alog(age) and 0.20 in Alog(Z), and 
merge adjacent bins when the noise level of the 
solutions indicates there is insufficient information 
content in the CMD to resolve the fine age bins. 

The CMD is divided into a regular grid of color- 
magnitude cells, and the expectation value of the 
number of stars in each cell for a SFR of 1 Mq 
yr~^ is calculated from the isochrones. There are 
several parameters that are taken as fixed con- 
stants during the solution. These include the 
IMF, fraction of binaries and mass ratio distri- 
bution function, and the distance mo dulus and 
redden ing. The adopted IMF is from IChabrier 



(|2003[ ). and the binary frac tion and mass ratios 
are pa ram etrized based on iDuquennov fc Mavor 
(|l99l[ ) and iMazeh et all (|l992l) . In this prescrip- 
tion, 35% of stars are single and the rest are bi- 
nary. The binaries are divided into "wide" and 
"close" binaries in a 3:1 ratio; the secondaries in 
the wide systems are drawn from the same IMF 
as the primaries, but in the close systems the sec- 
ondary masses are drawn from a flat IMF. The 
distances and reddeni ngs are initially co nstrained 
to the values given in Sabbi et al. I (l2009l) . but are 
varied if the resulting synthetic CMDs are mis- 
matched to the data. 

No age-metallicity relation is explicitly as- 
sumed, but a range of metallicities at each age 
is allowed, constrained by the color range of the 
data. In some cases, there is little leverage in 
the CMD to constrain the metallicity, so out- 
side information is used to choose the metallicity. 
For example, the V— I color of the upper MS is 
not strongly metallicity-dependent, so we adopt 



a metallicity of Z = 0.004 based on HII region 
and Cepheid metallicities for the youngest stellar 
populations in the SMC Bar. For stars with ages 
on the order of wl-7xl0* yr, this metallicity also 
gives a good match between the colors of the blue 
supergiant stars in SFHl and the Padova models. 
For older ages, ranges of metallicity at each age 
are allowed. 

A synthetic CMD is constructed by convolving 
the weighted, color-magnitude binned isochrones 
with the color and magnitude errors and incom- 
pleteness functions determined by artiflcial star 
tests. In the fit process, linear combinations of the 
individual synthetic CMDs are added to find the 
composite CMD that best matches the data. Be- 
cause many cells in the Hess diagram are empty 
or contain few stars, a maximum likelihood test 
based on the Poisson distribution must be used in- 
stead of a statistic (jCashlll979l) . A direct search 
of parameter space is infeasible because there may 
be dozens of age and metallicity bins in the so- 
lution; additionally, effects such as age-metallicity 
degeneracy can easily produce a large number of 
false, local maxima in the likelihood space. Be- 
cause of this we use a simulated annealing ap- 
proach, in which a simplex of initial guesses at 
the SFH is randomly perturbed and the changes 
are rejected with some probability if they worsen 
the fit. The star formation rates are transformed 
according to an arcsin function prior to fitting in 
order to prevent negative star formation rates from 
being considered. 

Errorbars on the SFH are calculated by test- 
ing each age-metallicity bin of the best-fit solu- 
tion in turn. The SFR in the bin under considera- 
tion is forced away from the optimal solution and 
the fitting procedure is redone, with the tested 
bin held fixed. The Itr error on the SFR of the 
bin is taken to be the value beyond which no flt 
can be found that is within Icr of the global best 
fit. Because the total number of stars is flxed, 
a deficiency of stars in one age bin can often be 
partially compensated by increasing the SFR in 
nearby bins; this means that errorbars are fairly 
large and the SFR in adjacent bins can be strongly 
anticorrelated. Because there are always unmod- 
elled populations (including Galactic foregrounds, 
background galaxies, and simple bad data) and 
there may be poorly- fit stellar sequences (e.g., the 
colors of RGB stars) , the overall fit quality as mea- 
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sured by the equivalent of a statistic is usually 
found to be quite poor; only the relative likeli- 
hoods are of any meaning. The most likely SFH 
to match the observed CMD is driven by the most 
populous cells in the CMD; because these are fre- 
quently near the faint end of the data where in- 
completeness is high, it is absolutely essential to 
have a very accurate model of the incompleteness 
to obtain robust results. 

The ACS data are of high quality, but incom- 
pleteness sets in at a relatively bright level. The 
large number of stars and low average error leads 
us to use magnitude bins 0.08 mag high and color 
bins 0.04 mag wide. In order to avoid systemat- 
ics due to the incompleteness, we restricted the 
comparison to magnitudes 15 < I < 22. Because 
there are significant uncertainties in modeling the 
colors of RGB and RC stars, we restricted our 
fits to the MS and SGB. We considered a set of 
isochrones binned by 0.10 in log(age) from 4 Myr 
< age < 13.5 Gyr (6.60 < log(age) < 10.13); below 
log(age) = 8.60 the decreasing number of bright 
tracer stars led us to double the size of the age 
bins. We considered metallicities of Z = 0.00015, 
0.0004, 0.0006, 0.001, 0.0015, 0.0024, and 0.004. 

We begin with the distance moduli and red- 
denings from ISabbi et al. I dlooi), but found that 
these needed to be modified in order to obtain 
the best match to the data with the Padova 
isochrones. For SFHl we found the best distance 
modulus (m-M)o = 18.83 and reddening E(B-V) 
= 0.08, while in SFH4 we used 18.85 and 0.12. In 
both fields, the values are within la of the initial 
guesses. In SFHl we found that the model upper 
MS was significantly narrower and bluer than the 
data with the adopted reddening value, so we were 
forced to introduce an extra reddening component. 
We adopted the simple expedient that all popula- 
tions younger than log(age) — 7.60 were subject 
to double the mean reddening of the field, while 
stars with 7.60 < log(age) < 8.00 were assigned a 
reddening halfway between the younger and older 
stars. Differential reddening in th e central SMC 



has p reviously been observed by iZaritskv et al 



(|2002[ ). who used multicolor photometry to derive 
line of sight extinction corrections to individual 
bright SMC stars. They found a mean Ay — 0.18 
mag for stars with 5500 K < Te// < 6500 K, and 
Ay = 0.46 mag for stars with Te// > 12000 K; 
these numbers are in reasonable agreement with 



the va lues adopted here for SFHl. IZaritskv et al 
(j2002n discuss the physical reality of their derived 
reddening distributions and conclude that small 
scale reddening variations may reasonably be at- 
tributed to residual gas and dust in star-forming 
regions. The vast majority of old, cool stars will 
not be seen through star-forming regions, and so 
the small highly- reddened tail of the old popula- 
tion will not strongly influence the mean redden- 
ing. Because the amount of differential reddening 
has high spatial frequency, the numerical agree- 
ment (as in SFHl) or lack thereof (as in S FH4) 



with the conclusions of lZaritskv et al. I (|2002h must 



be considered fortuitous. 

Cole's best-fit solutions are presented here as 
SFHl-C and SFH4-C respectively for the two 
fields. 

4. Results: SFH of two fields in the SMC 
Bar 

The properties of the three SFH solutions for 
each field are summarized in Table [TJ While 
the two procedures use the same photometry and 
isochrones, they differ in a number of important 
respects. The application of two completely in- 
dependent codes to derive the SFH allows us to 
identify the most model-independent features of 
the SFH, gives insight into the areas of the CMD 
that are driving the fitting procedures to the re- 
sults they give, and gives an estimate of the sys- 
tematic uncertainties in the fits. 

4.1. SFHl field 

The SFHl-A (restricted metallicity) solution 
which gives the lowest (2.6) reduced shown 
in Fig. |4] (error bars represent 1-cr uncertainty), 
while the corresponding synthetic CMD is shown 
in the middle panel of Fig. [5l The left panel 
of the same figure shows the observational CMD 
for a direct comparison. Our best distance mod- 
ulus and reddening are (m — M)q — 18.77 and 
E{B — V) = 0.11, respectively. According to this 
solution, the SFR has been extremely low during 
the earliest 6-7 Gyr, with a strong and rapid in- 
crease around 5 Gyr ago. From then on, the star 
formation activity has been typically gasping, with 
ups and downs with respect to an average almost 
constant rate. 

In the middle row of Figure [6] we plot the lu- 
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Table 1: Summary ol SFH Solution Parameters" 



Solution 


Method 


(ni-M)o 


E(B-V) 


IMF 


CMD binning 


Metallicities 






(mag) 


(mag) 




(color X mag) 


(ZxlO^) 


Field SFHl: 














SFHl-A 


Bologna 


18.77 


0.11 


Kroupa (2001) 


variable'' 


1, 2, 4 


SFHl-B 


Bologna 










0.4, 1, 2, 4, 8 


SFHl-C 


Cole 


18.83 


0.08"^ 


Chabrier (2003) 


0.04x0.08 


0.15, 0.4, 0.6, 1.0, 1.5, 2.4, 4.0 


Field SFH4: 














SFH4-A 


Bologna 


18.80 


0.11 


Kroupa (2001) 


variable'' 


1, 2, 4 


SFH4-B 


Bologna 










0.4, 1, 2, 4, 8 


SFH4-C 


Cole 


18.85 


0.12 


Chabrier (2003) 


0.04x0.08 


0.15, 0.4, 0.6, 1.0, 1.5, 2.4, 4.0 



"All models based on the Padova isochrone set; see text for details. 
''See Figure[3] 

^Differential reddening assumed for stars younger than 100 Myr. 



minosity function (LF) of blue and red stars (first 
and second column from the left, respectively) in 
the observed CMD and in the one resulting from 
SFHl-A. The top and bottom panels show for 
comparison the result when the SFH is searched 
using IMF exponents s = 2.0 and 2.5, respectively, 
for stars above 0.6 Mq. 

Despite a somewhat high reduced 2.6, this 
model is successful in describing the main features 
of the data. 

However when the results are examined in 
detail a few noticeable issues can be identified. 
First of all the synthetic upper MS is sharper 
tha n the obser v ationa l one. As already noted 
by ISabbi et all (|2009h . such a broadening may 
indicate an additional absorption. One explana- 
tion could be that the SFHl field is in the main 
body of the SMC and the reddening material may 
be patchy and cover not uniformly all lines of 
sight. Indeed, this is what is obs e rved in the 
extensive analysis of iHaschke et alJ (|201ll ) , who 
suggests for SFHl a broader reddening distribu- 
tioi|l( cr^(y_/) w 0.12 mag) with respect to SFH4 
{cte{v-i) ~ 0-09 mag). Yet other possibilities are 
that a fraction of massive MS stars is affected by 
rotation, which can lead to widened upper M Ss in 

or a 



young clusters (see e.g. Grebel et al. 1996[ ). 



different SMC depth along SFHl and SFH4's line 
of sights. 

Our model slightly but systematically underes- 
timates the number of blue massive stars brighter 
than V ^ 19 by 20—30 percent (see the middle 
panel of the first column in Fig. |6]). A fiatter 



•^Reddening values from the red clump method are available 
through the German Astrophysical Virtual Observatory in- 
terface at |http://dc.zah.uni-heidelberg.de/mcx| 



IMF could mitigate such discrepancy (see the top 
panel of the first column in Figure ^ , but the cor- 
responding model would underestimate the num- 
ber of low mass stars (see the discrepancy between 
V = 2A and = 26 in the top panel of the second 
column in Fig. |6]). 

As far as the evolved stars are concerned, the 
synthetic CMD matches well both the SGB magni- 
tude spread (that is a signature of the age spread) 
and the RGB color dispersion below the RC lumi- 
nosity. However, there are several differences as 
well: 1) the predicted RC morphology is irregu- 
larly shaped, while the observational RC is smooth 
and rather elliptical. This difference may be par- 
tially explained by the coarse metallicity resolu- 
tion of our model, but also by a small amount of 
differential reddening; 2) the synthetic CMD over- 
estimates the number of RC stars by about a fac- 
tor of two; a steeper IMF could slightly mitigate 
this mismatch (see the bottom panel of the second 
column in Fig. [6]) but at the expenses of the up- 
per MS fit; 3) the RGB stars brighter than the RC 
are too blue (even if the predicted counts match 
exactly the observed ones) . This suggests that our 
models are too metal poor or, more likely, that the 
adopted color transformations systematically fail 
near the RGB tip. 

We investigate the effect of metallicity on the 
SFH found via the Bologna procedure by adding 
the Padova tracks with Z = 0.008 and Z = 0.0004 
(solution SFHl-B). The top-left panel of Fig. [7] 
shows the resulting SFH, while the other panels 
of the same figure present the mass fraction con- 
tributed by each metallicity. 

Overafi, the solutions SFHl-A and SFHl-B are 
rather similar (see Fig. |8l), both being character- 
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Fig. 4. — SFH for the SFHl field, obtained with a restricted metaUicity range (solution SFHl-A). The 
star-formation rate is given both in Mq yr~^ (left ordinate) and per unit area (Mq yr^^ kpc^^, in the right 
ordinate). 



ized by a long quiet period at the oldest epochs. 
However, solution SFHl-B shows some different 
features. First, the early activity is slightly higher, 
a direct consequence of the lower initial metaUic- 
ity; then, the activity between 0.5 and 3 Gyr ago 
is smoother than in SFHl-A. 

We find with additional test models that the 
epoch of the first peak of star formation activity 
is progressively earlier for decreasing values of the 
adopted lowest metaUicity. We are thus confident 
that our conclusion of a long almost quiescent pe- 
riod earlier than 5 Gyr ago is robust, although we 
do see and predict some stars as old as 10-12 Gyr. 

Concerning the recovered chemical history, we 
find that the largest stellar mass fraction is pro- 
duced at Z = 0.001, while the contribution from 
Z = 0.004 is only a few percent. Generally, the 
metaUicity increases with time. 

From the point of view of the fit quality, the 
new solution improves the CMD match, yielding a 
reduced oi about 2.3. The differences are vis- 
ible both in the synthetic CMD of Fig. [S] (right 
panel) and in the LF (third and fourth columns) 
of Fig. [6] Model SFHl-B reproduces the upper 
MS star counts and morphology better than model 



SFHl-A (compare the middle panel of the third 
column with the middle panel of the first column 
in Figini), but some problems stUl affect the sim- 
ulation: 1) the synthetic CMD shows a hint of a 
red HB, which is not observed in the data; 2) the 
number of BL stars is still under-predicted; 3) the 
RGB above the RC is too blue; 4) the predicted 
RC is still overpopulated and irregularly shaped. 

Given these results, especially the improvement 
along the upper MS, model SFHl-B is globally 
better than SFHl-A, even if neither of them is 
fully satisfactory. 

The result of applying Cole's CMD-fitting pro- 
cedure to SFHl is shown in Figure |9l The 
left panel shows the data, binned 0.04x0.08 in 
(uiy—uii, m/), while the middle and the right 
panels show the best-fit synthetic CMD (SFHl-C) 
and map of residuals, respectively. The fit proce- 
dure has matched the LF and mean color of the 
stellar sequences well in general. The data and 
the model can be readily distinguished from one 
another because the data is just one instance of 
a random draw from the parent population and 
contains unmodelled noise, while the model repre- 
sents the best guess at the pure parent population 
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Fig. 6. — SFHl: Predicted (red dot-dashed histogranis) and observational (black solid histograms) LPs for 
stars bluer (first and third columns) and redder (second and fourth columns) than V — I = 0.6 for the model 
solutions SPHl-A (first two columns) and SPHl-B (third and fourth colums). Models in the top, middle 
and bottom panels are computed with IMP exponents above 0.6 Mq s = 2.0, 2.3, 2.5, respectively. 



and is therefore more smoothly distributed. 

Among the well-matched features are the color 
and vertical extent of the RC and its upward ex- 
tension to I «17, the enhancement in the SGB 
at I w20.8, and the nearly-vertical finger of stars 
at V— I wO.5 corresponding to an intermediate-age 
MSTO. Notable areas of mismatch include the fail- 
ure to reproduce the width of the MS (most ob- 
vious for I<21), and the factor of 2 overproduc- 
tion of RC stars (I wl8.5, V-I wl). There is a 
sparsely populated red HB in the model that is 
not apparent in the data; this may be related to 
the general factor of 2 overproduction of low-mass 
core helium-burning stars, to poor constraints on 
the detailed element abundances at ancient times. 



and/or to gaps in modeling the physics of HB en- 
velopes. 

The SPH over the period from 0.5 to 13.5 Gyr 
ago for SPHl-C is given in Figure [TUl The SPH is 
characterized by a very low star formation rate for 
several billion years before a rapid increase about 
5 Gyr ago. This is the event that produces the 
features in the MS and SGB of the SFHl CMD. 
The SPH remains at a similar elevated level after 
the 5 Gyr event, with some fluctuations. 

OveraU, the SPH of the SFHl field, based on 
the synthetic CMDs from the Bologna and Cole 
procedures, is characterized by the following fea- 
tures: 
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Fig. 7. — SFH for the SFHl field resulting from models which include also the metallicities Z ~ 0.008 
and Z = 0.0004 (solution SFHl-B). The top-left panel shows the total SFH while the others display the 
contribution from each of the labeled metallicities. 



• The first 6-8 Gyr were rather quiescent and 
only a small fraction of the stars in this 
field was formed during these old epochs. 
The inferred average rate of star forma- 
tion for ages older than w5 Gyr is only w3 
xlO"^ Mq yr^i kpc"^ 



About 5-6 Gyr ago SFHl experienced a re- 
markable enhancement in the stellar produc- 
tion: over wl Gyr the activity ramped up to 
about 2.2x 10"^ M© yr'^ kpc"^. The age 
and magnitude of this enhancement is ro- 
bust to choices of IMF, CMD gridding, red- 
dening, assumed metallicity, and details of 
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Fig. 8. — Recovered SFHs for SFHl using the Bologna procedure: SFHl-A vs. SFHl-B. The latter is a 
better match to the data. 



the fitting procedure. 

• The average SFH has not dropped signif- 
icantly from its peak at 5 Gyr ago, but 
the degree of burstiness in the solutions is 
model-dependent. SFHl-B shows a rela- 
tively smooth recent history, while SFHl-C 
shows a factor of two drop in SFR between 
3-4 Gyr ago with a subsequent recovery and 
SFHl-A shows repeated gasps from 0-2 Gyr 
ago. The reasons for this range of behavior 
are considered in section l4?3l below. 

It is worthwhile to quantify the stellar mass 
produced in each age interval (see the cumulative 
mass function in Fig. [TT|) . According to solution 
SFHl-B, the SFHl field assembled a small fraction 
(11%) of its stellar mass before 6 Gyr ago, while 
a significant fraction (about 40%) was assembled 
only in the last 2 Gyr. For comparison, SFHl-C 
predicts that 17% of the stars were formed prior 
to 6 Gyr, and 33% over the past 2 Gyr. 



4.2. SFH4 field 

The history of the SFH4 field was determined 
in similar fashion to the SFHl field, using three 
sets of simulations (see Table [T]). The recovered 
SFH for case SFH4-A is shown in Fig. [H 

As in SFHl, the star formation proceeded at 
a low level until 5-6 Gyr ago, when it rose to al- 
most the same amplitude of the first peak in SFHl. 
Afterward, the actvity started a slow but steady 
decline until now and only in the last 50 Myr it 
reached higher levels (mimicking in this respect 
the behavior of SFHl-A). Concerning distance and 
reddening, the recovered values ((771— M)o = 18.80 
and E{B — V) = 0.11) are not significantly dif- 
ferent from those obtained for SFHl. The cor- 
responding reduced about 1.8, is significantly 
lower than the best value obtained for SFHl. In- 
deed, a visual inspection of the synthetic CMD 
(middle panel of Fig. I13|) confirms a very good 
agreement with the observational counterpart: our 
best model can reproduce the position and the 
morphology of the MS, the SGB and the RGB 
(above and below the RC). The only minor mis- 
matches are in the BL region, which is underpop- 
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Fig. 9. — Left panel: observed Hess diagram for the SFHl field of the SMC Bar, created by binning the data 
by 0.04 in color and 0.08 in magnitude; middle panel: synthetic CMD from the SFHl-C model (see text for 
details); right panel: difference Hess diagram between the data and SFHl-C solution. 



ulated in our model, and in the SGB luminosity 
distribution, which is slightly more discontinuous 
than in the data. 

As for SFHl, we expanded the metallicity range 
to include the values Z = 0.008 and Z = 0.0004 
and re-derived the SFH. The essential features of 
the new solution (SFH4-B; see Fig. IM)) are not 
changed (see Figure [15] for a comparison SFH4- 
A vs SFH4-B): there is still a long quiet initial 
period, followed by a rapid enhancement of the 
activity and a subsequent decline. According to 
the solution SFH4-B, most of the stars ever formed 
in SFH4 have metallicities around Z = 0.001. 

In terms of fit quality, solution SFH4-B im- 
proves the match to the data significantly, leading 
to a reduced of about 1.3. The corresponding 
synthetic CMD is excellent fFigflSl right panel) 
and the only detectable discrepancy is around the 
BL phase, still slightly underpopulated by the 
model. Unlike in SFHl, the number of RC stars is 
well matched (see the LF in the right panel of Fig. 
ire]) . The fit to SFH4 using Cole's method shows 
similar results, with a better fit to the upper MS 
than in SFHl because the observed sequence is 



narrower in color, and an overproduction of RC 
stars. The observed Hess diagram, synthetic di- 
agram from solution SFH4-C, and the difference 
between the two are shown in Figure [TT] 

SFH4-C is broadly similar to SFHl-C, but the 
mean SFR is lower at all ages, commensurate 
with the smaller number of stars in the field (Fig- 
ure ITSl) . It is notable that the onset of significant 
star formation occurs simultaneously in the two 
fields, to within the ±1 Gyr precision of our data. 
However, the SFH4 field shows a decrease in SFR 
relative to SFHl over the last «1.5 Gyr (which is 
obvious from a comparison of the RC morphology 
and upper MS LF of the two fields). 

According to the solution SFH4-B (see the cu- 
mulative mass function of Figure [T9| about 60% of 
the stellar mass has been produced between 6 Gyr 
and 2 Gyr ago, which is higher than the fraction of 
mass (48%) produced in SFHl in the same period. 
On the other hand, 16% of the total mass was al- 
ready in place before 6 Gyr ago, which is compara- 
ble with the production in SFHl (11%). The for- 
mer difference is compensated in the last two Gyr, 
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Fig. 10. — Long term history of the SFHl field from sohition SFHl-C. Significant star formation commenced 
«5 Gyr ago and the rate has remained high, with some fluctuation, ever since. 



when SFHl has turned into stars about twice the 
mass fraction of SFH4. While in the SFHl field 
the Bologna solution produced a younger result 
than the Cole solution, in SFH4 the situation is 
reversed. The solution SFH4-C predicts that only 
7% of the stellar mass was in place before 6 Gyr 
ago, while 29% was astrated in the past 2 Gyr. 

4.3. Comparison of the two methods 

We have derived the SFH of the field stars in 
the SMC Bar using two different synthetic CMD 
techniques. Figure [20] shows the results plotted 
together on the same scale to illustrate the simi- 
larities and differences. For SFHl both SFHs indi- 
cate a major activity in the last 5 Gyr, an overall 
plateau (on average) since then on, and a spike 
in the last 50 Myr. The major differences con- 
cern: 1) The exact epochs of the peaks between 
250 Myr and 4 Gyr ago (around 1.5 Gyr ago in 
Bologna's solution and 2.5 Gyr ago in Cole's so- 
lution); 2) The rate in the period 8-13 Gyr ago, 
which is slightly stronger in Cole's solution. 

The agreement for the field SFH4 is good as 
well: both models predict a star formation onset 



between 4 and 5 Gyr ago, followed by a slow de- 
cline and a very recent burst. The major difference 
here is the rate of decline in SFR after the 4-5 Gyr 
enhancement, with Cole's solution showing a con- 
stant activity for about 2.5 Gyr and Bologna's so- 
lution showing a faster decline followed by a mild 
enhancement around 1.5 Gyr ago. 

In both fields, the only ages in which the star 
formation rates differ by significantly more than 
the formal errorbars on the solutions are in the 
period 1.5-3 Gyr ago, with Cole's solution show- 
ing consistently higher SFR. Part of this differ- 
ence stems from the slightly higher SFR derived 
by Cole's method at all ages, owing to the different 
mean stellar masses resulting from the different 
assumed IMFs and binary fractions. The effect 
may be exacerbated in the 1.5-3 Gyr age range 
because of the interplay between age and the dif- 
ferent metallicity values considered, and the fact 
that Cole's procedure does not consider the RC 
in calculating the best-fit SFH. Consideration of 
a very large number of metallicity values tends to 
produce smoother SFHs tha n those derived u sing 
only a few metallicities (e.g.. ICole et al. I I2OO9I) . 
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Fig. 11. — Cumulative mass function in SFHl according to the solutions SFHl-B and SFHl-C. 



Concerning the distance modulus, Cole's best 
values (18.83 for SFHl and 18.85 for SFH4) are 
slightly higher than Bologna's (18.77 and 18.80). 
In SFHl these offsets may be due to the higher 
reddening required by the Bologna method, but 
in SFH4 the converse is true. Yet the distance 
modulus offsets are similar. The offsets are quite 
small and can be considered to be within the 
noise, but a possible reason could be the fact 
that Cole's method ignores the RGB and RC, 
so the distance is essentially a MS-fitting dis- 
tance, while the Bologna method includes RGB 
and RC information. More interestingly, both dis- 
tance estimates are shorter than recent determi- 
nations based on RR -Lyrae ((to — Af)o — 18.90, 
iKapakos et al. 2011) and eclipsin g binaries ((m — 
M)o = 19.11. iNorth et al.ll2010l ). but still com- 
patible with the average distance of star clusters 
(around 18.87 for lGlatt et"aLll2008bl and between 



18.71 and 18.82 for lCrowl et al.ll2001 ). This vari- 



ance may indicate a different sensitivity to redden- 
ing, wln£hj£ju£hly_v^ (see 
e.g. lZaritskv et al.l[20ollHaschke et al.ll201ll ). or a 
lin e of sight depth effect (up to 4.9 kpc according 
to Subramanian &: Subramaniam 2009i between 



10 and 17 kpc according to iGlatt et al.l l2008b[ ) 
which may depend on the direction. 

The age-metallicity relations inferred by both 
methods are consistent, implying a metallicity 
slowly increasing with time and the bulk of the 
stars with metallicities between 0.001 < Z < 0.002. 
The models derived with the Bologna procedure 
consistently find better fits when a metallicity of 
Z = 0.008 is included for the younger ages, while 
such a population is not required in Cole's models 
(see Fig. [21]). However, Cole adopted a higher dif- 
ferential reddening in his modeling of SFHl, which 
has the higher fraction of young stars; see Table[T] 

4.4. Comparison Between the Two Fields 

Spatially Resolved SFH of the SMC Bar 

Direct comparison of the two fields shows that 
colors and magnitudes of the prominent stellar se- 
quences are similar, and the main differences con- 
cern the younger stellar populations. In SFHl the 
upper MS, BL, and red supergiant region are more 
populous, indicating a higher SFR at ages of a few 
hundred Myr or less. A challenge for both SFH- 
fitting methods was to reproduce the breadth of 
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Fig. 12. — Best recovered SFH with restricted metallicity range (solution SFH4-A) for the field SFH4. 



the upper MS in SFHl, as well as its RC, which 
was overproduced by both procedures. Integrat- 
ing up the predicted stellar masses from the SFH- 
fitting procedures we find that «4xl0^ Mq stars 
formed over the lifetime of SFHl, and 2x10^ Mq 
in SFH4. 

The general similarity of the two CMDs indi- 
cates that their features are representative of the 
general properties of the central SMC over most 
of its history. Small localized bursts/gasps of star 
formation may be smoothed away by the drift- 
ing of stars away from their birthplace except for 
the youngest ages. The dramatic increase in SFR 
at 5 Gyr is a strong feature in both fields and 
suggests a global change in environment to one 
that strongly favored star formation throughout 
the SMC at this time. Backward integration of 
the Magellanic Clouds orbital paths around each 
other and the Milky Way over this length of time 
are quite uncertain, and it is therefore difficult to 
say if this increase is associated with a tidal inter- 
action or not. 

More recent dynamical signatures should be 
easier to locate in time, and thu s to correlate 
with global star formation events. iDiaz fc Bekki 



( 201l[ ) presented orbital calculations for the sys- 
tem MW/LMC/SMC, finding evidence for a close 
encount er about 1.5- 2.0 Gyr ago. Analyzing SMC 
clusters JPiattil (|201l[ ) also suggested two enhanced 
formation processes that peaked at 2 and 5-6 Gyr 
ago. A comparison with our solutions suggests 
that the former episode might be associated with 
the Bologna (1.5 Gyr) and Cole (2.5 Gyr) inter- 
mediate age peak, while the latter coincides with 
our strong rise 5 Gyr ago. 

The inset panels in Fig. [20] show the recent 
1 Gyr history of both fields according to Bologna 
and Cole methods. In the SFHl field the mean 
specific SFR over the past 1 Gyr is nearly as high 
as the Gyr-averaged peak rate during the 5 Gyr 
event, but in SFH4 the recent SFR is reduced by a 
factor of ~3 from the intermediate-age peak. How- 
ever, in both fields the recent average specific SFR 
is still several times higher than the rate prior to 
the major episode 5 Gyr ago. The increase in the 
very youngest age bin may be an artifact of the fit- 
ting process, because the only stars younger than 
«50 Myr in our CMDs are still on the MS - all 
of the evolved stars are above our bright limit- 
and so the code is free to vary the SFR at young 
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Fig. 13. — Synthetic CMD for the solution SFH4-A (middle panel) and SFH4-B (right panel) compared with 
the observational CMD (left panel). 



ages somewhat arbitrarily to fit the residuals left 
after older populations are constrained. However, 
in literature there are also reports indicating very 
recent bursts of activi t y; see e.g. the analysis by 
Indu fc Subramaniaml ( 201l[ ) who found peaks at 
0-10 Myr and 50-60 Myr ago. Stars formed in 
these events, however, may not have had time to 
diffuse through the SMC, and thus may not be 
present in our two small HST fields. 

There are also clear differences between the two 
fields, with SFH4 showing a lower overall SFR den- 
sity averaging wSxlO^^ M© yr~^ kpc~^ and SFHl 
forming stars at the average rate of 2.8xl0~^ 
Mq yr~^ kpc^'^ with an apparently higher de- 
gree of burstiness. The features in the CMD that 
indicate high SFR at the "burst" ages are the 
large number of MS stars at I wl7-17.5 (100- 



150 Myr) and the steepening of the MS LF at 
around I wl9.5-20, with the resulting large num- 
ber of red core helium-burning stars (vertical RC; 
400-800 Myr). In SFHl the average recent SFR 
is comparable to the peak of the longterm average 
SFH since the 5 Gyr episode, but in SFH4 the rate 
has dropped by a factor of from the peak. If 
it is proposed that a t idal interaction be tween the 
SMC and LMC (e.g. iBesla et al.ll2012l ) or Milky 
Way resulted in a period of increased SFR at 100- 
200 Myr as seen in SFHl, then it remains to be 
explained why this event is not seen in the SFH4 
field. 

4-4-^- Chemical Evolution 

Because a range of metallicities is considered 
at most ages, we can plot the resulting age- 
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Fig. 14. — Best SFH (solution SFH4-B) resulting from models which include also the metallicities Z = 0.008 
and Z = 0.0004. The top-left panel shows the total SFH while the others display the contribution from each 
of the labeled metallicities. 
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Fig. 15. — Comparison between the solutions SFH4-A and SFH4-B. 



metallicity relation for the best-fit solutions. This 
is shown for the Cole solutions (which use a finer 
grid of metallicities, see Table [IJ in Figure [HJ 
which shows that the two fields have experi- 
enced similar longterm chemical evolution his- 
tories. This is not a surprising result because 
over time periods longer than ^10^ yr the entire 
galaxy should be mixed, but it is a confirmation 
that the fit procedure is giving consistent results 
despite the differences between the two fields. For 
the youngest populations, only a single metal- 
licity (Z=0.004) has been used, so the errorbars 
simply reflect the spacing of the metallicity grid. 
Where multiple metallicities were used, the mass- 
averaged mean metallicity is given, and the error 
bar shows the range of metallicities within which 
90% of the stars in each age bin are expected 
to fall. The mean metallicity is quite similar to 
that obtained from the Bologna method results, 
with the exceptions that Z=0.008 populations are 
not present in large numbers in the Cole solution, 
and the Bologna solution for SFH4 appears to in- 
crease its metallicity more quickly for ages older 
than «6 Gyr. Given the uncertainties in RGB 
model colors and the differences in distance and 



reddening between the two fits, it is not clear these 
differences are significant. 

The field star AMR sh own in Fig. [5T] is con- 
sistent with IPiattil 's (|2012l ) AMR, whose AMR is 
obtained from CCD Washington CTl photome- 
try. Moreover the metallicity values for ages older 
than 1 Gyr match well the spectroscopic field star 
measureme nts of red giants in the central SMC 
reported bv lParisi et al. I (I2OIOI) . 

Compared to AMRs derived from clusters, our 

result is in go od agreemen t with Da Costa fc Hatzidimitrioul 



(119981) and iGlatt et al.l (l2008bD. w hile it dif 



fers from the clus t er AM R of iPiattil (|2011l) and 
Parisi et al.l (|2009f ) . IPiattil 's (|201l[) metalhcitv dis- 



persion is much higher than our formal errors at 
any age, suggesting that SMC clusters may have 
originated fr or n a le ss mixed interstellar medium. 
Parisi et al.l 's ( 2009( ) AMR rises faster than ours 



prior to 10 Gyr, suggesting that SMC clusters 
may have experienced an independent chemical 
evolution history. 

It is worth to remind that we have derived the 
best AMRs from two small regions of the SMC 
and more fields are needed to confirm our results. 
A detailed comparison of our six AMRs (sampling 
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Fig. 16. — Predicted LF (red histograms) for the solution SFH4-B vs the observed (black histograms) LF. 
The left-hand and the right-hand histograms are for stars bluer and redder than V ~- I — 0.6 respectively. 



regions of the Bar, Wing and Halo) with spectro- 
scopic metallicity measurements will be presented 
in a future paper. 

5. Comparison with other studies 



Harris fc Zaritskvl (|2004l ) (hereafter HZ04) were 



the first to apply the synthetic CMD method to 
the derivation of the SMC SFH. They used wide- 
field UBVI photom etry of an area of 4° x 4.5° 
(MCPS catalog. Zaritskv et al.l[2002[ ). suggesting 
that 50% of SMC stars formed earlier than 8.4 
Gyr ago, and that the period between 3 and 8.4 
Gyr ago experienced very little, if any, star for- 
mation activity. Figure [22] shows HZ04's results 
for their regions including our SFHl and SFH4 
compared with our derivations. Both HZ04's so- 
lutions show a significantly higher stellar produc- 
tion prior to 8.4 Gyr ago. At intermediate ages 
HZ04's solution for SFHl is characterized by a 
clear and long lull between 2 Gyr and 8.4 Gyr 
ago, which is in striking contrast with our SFHs, 
whereas HZ04's SFH for SFH4 is stronger and in 
satisfactory agreement with our predictions. In 
the last 2 Gyr both HZ04's star formation rates 
are much stronger than our findings. 

One should note that the HZ04 SFHs are de- 
rived from a much larger field of view (about 
16 times wider) and their photometry did not 



reach the oldest MSTO. This last point is partic- 
ularly relevant, since studies that do reach it indi- 
cate, instead, that, although present, stars older 
than 8 Gyr do not doni i nate the SMC popula- 
tion toolDhin et al.ll200l t iMcCumber et al.ll200£ " 



Noel et all 12007. 2009; C hiosi fc Vallenar: 



a]l2005l: 
fl2Q07|). 



Further support to this is also provided by the rel- 
atively low number of RR Lyrae stars detected in 



the S MC compared to the LMC (jSoszvhski et al 



2OIOI) . 



Dolphin et al 



(|2001[ ) analyzed an external re- 
gion close to the globular cluster NGC 121, and 
concluded that stars in the outskirts of the SMC 
formed during a broadly peaked episode of star 
formation, with the largest rate between 5 and 8 
Gyr ago. As expected, this outer field is signifi- 
cantly older than ours, although the contribution 
from the age bins older than 10 Gyr is again a 
minor fraction of the total mass produced. 

studied an external 
found an increas- 
and then over the 



McCumber et al 



(120051) 

Wing, and 



field in the SMC 
ing rate from 12 to 4 Gyr ago, 
past 1.5 Gyr, with a significantly quieter phase 
between 4 and 1.7 Gyr ago. This analysis was 
not conducted using a statistical significance test- 
ing and their favored solution was the one (among 
three ansatz models) which best matched the star 
counts in five strategic CMD regions. Although 
this makes the comparison with our finding diffi- 
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Fig. 17. — Left panel: observed Hess diagram for the SFH4 field; middle panel: synthetic Hess diagram from 
solution SFH4-C; right panel: difference between the two. 



cult, two considerations can be made: 1) also in 
this case the stellar production prior to 10 Gyr ago 
is low, in agreement with our result; 2) apart from 
a recent burst, the average activity of their SFH is 
earlier than in our fields, as expected for a Wing 
region, which represents a population with inter- 
mediate age between the Bar and the outskirts. 

Noel et al.l ([2009) 



Noel et al. 



(|2007n and llNoel et al.l ([^UUyj pro- 
duced one of the most extensive and accurate 
photom etric campaign i n the SMC outskirts (see 



also Nid ever et al.l 120111 ). Their 12 fields are suf- 



ficiently deep to reach the old MSTO and none 
of them showed a clear extended HB, representa- 
tive of a very old and m etal poor stellar popu- 
lation. Noel et al. I (l2009l) used a population syn- 



thesis technique and recovered a great variety of 
SFHs, with only the outer fields characterized by 
a strong old and intermediate age activity. Al- 
though the spatial distribution of these fields, 
mainly enclosed in the periphery of the SMC, 
makes the comparis on with our results less direct, 
we note that among iNoel et all 's (|2009h solutions, 
their innermost field (qj0112) most closely resem- 
bles our history of the Bar. Prior to 8 Gyr ago, the 



typical SFR density in their fields seems compara- 
ble with ours, suggesting that the SMC periphery 
and Bar share a common old population but that 
the initiation of significant star formation at 5-6 
Gyr affected the Bar far more than the external 
regions of the SMC, marking the birth of the Bar 
as a distinct feature. 

A more direct c o mpar ison can be made with 
Chiosi fc Vallenari (l2007h . whose analysis was 



focused on three deep HST/ACS fields located 
in the Bar. In this case, our solutions are in 
good quantitative an d qua litative agreement. All 
Chiosi fc Vallenaril 's (|2007t ) SFHs show both an 



unambiguous rise between 7 and 5 Gyr ago (with 
the precise epoch varying depending on the spe- 
cific direction) and a negligible earlier activity. 
However, despite the low rate, their and our early 
activities are not zero, offering a natural expla- 
nation to the relatively low number of RR-Lyrae 
stars found in the SMC. Moreover, no two solu- 
tions are alike, suggesting that the star formation 
rate at young and intermediate ages was strongly 
variable across the SMC Bar. 
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Fig. 18.— Long term history of the SFH4 field from solution SFH4-C. Similar to SFHl, the field "switched 
on" «5 Gyr ago. The SFR has been declining over the past «1.5 Gyr in this locale. 



Our results are rather similar to the SFHs of 
other two Local Group irregulars such as IC 1613 



(ISkillman et al.l l2003r ) and Leo A (jCole et al 



20071 ): both galaxies show an initial quiet pe- 
riod followed by a prompt rise in the SF activity 
about 5-7 Gyr ago, with a peak between 3 and 6 
Gyr and between 2 and 3 Gyr ago, respectively. 
In both cases the recent activity (< 2 Gyr ago) 
is more similar to that in SFH4 than in SFHl, 
presumably because the IC 1613 explored field is 
in the outskirt and the Leo A field is large enough 
to include central as well as peripheral regions. 

Finally, comparing our SFH with the solutions 
found in the literature for the LMC (see e.g 
Harris fc Zaritskv 20091 and ISmecker-Hane et al 



20021 ) we find an overall similarity. Apart from 
the very early SF enhancement predicted by 
both Harris et al. and Smecker-Hane et al. 
that has no counterpart in our SMC solutions, 
our results for SFH l share key properties with 
Smecker-Hane et all 's (|2002[ l findings for the LMC 
Bar: 1) a lull over the period 6-10 Gyr ago; 2) a 
steep increase at 5-6 Gyr; 3) a gasping regime over 
the last 5 Gyr. 



In spite of th ese similarities, the SMC and LMC 
AMR s (see e.g. Carrera et al.ll2008 . Harris fc Zaritskv 
2009t ) remain rather different. While the LMC 
metallicity has increased faster both prior to 
10 Gyr ago and over the last 4 Gyrs (reaching 

[Fe/i?] 0.7 and [i^e/iJ] = -0.2, respectively), 

between these periods it has progressed at much 
slower rate. 

Moreover, it is well known that the LMC clus- 
ter distribution shows a long gap between 3 and 
12 Gyr. If this gap is a consequence of the qui- 
escent period between 6 and 10 Gyr ago (see e.g. 
[Harris fc Zaritskv 2009f ). then it is intriguing that 
a similarly long quiescent period has no counter- 
part in the SMC cluster distribution (which has 
been steadily increasing with time). 

6. Summary and conclusions 

This paper is the first of a series devoted to 
quantitative reconstruction of the SMC SFH from 
the deep HST/A CS observations presented in 



Sabbietal 



( 20091 ) ■ Here we have explored the 
directions SFHl and SFH4, both located in the 
SMC's Bar, by comparing the observational CMDs 
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Fig. 19. — Cumulative mass function in SFH4 according to the solutions SFH4-B and SFH4-C. 



with a library of model CMDs incorporating pho- 
tometric uncertainties and i ncompleteness as esti- 
mated bv lSabb i et al.l I (l2009l) . To provide a robust 
characterization of the SFH, the choice of the best 
model CMD was independently conducted with 
two objective sta tistical method s, namely Cole's 
(ICole et al.ll2007t ) and Bologna's (|Cignoni fc Tosi 
20101 ) procedures. 

Our best simulated CMDs exhibit an overall 
good agreement with observational CMDs. The 
star counts along the MS and the SGB morphol- 
ogy are generally well reproduced, indicating that 
our recovered SFHs and metallicity distributions 
are reasonable. However, while SFH4 CMD is well 
fitted, there are some difficulties to reproduce the 
exact morphology of SFHl's CMD, especially the 
upper-MS spread and RC counts, which are under- 
estimated and overestimated respectively. Con- 
cerning the resulting SFHs, a good consistency 
is found between the two methods in both fields. 
The only significant difference is the stronger rate 
suggested by Cole's SFH4 solution between 1.5 
and 3 Gyr ago. 

The combination of synthetic CMDs which 
most resembles the observations suggests the fol- 



lowing picture. At early times, both fields experi- 
enced a long quiescent phase characterized by low 
SFRs, followed by a rapid SF increase around 5-6 
Gyr ago. Since then, the mode of star formation 
has been somewhat different in the two fields. In 
SFHl the star formation was gasping and reason- 
ably high up to today. In SFH4 it was smoother 
and slowly declining. 

To account for these differences and similarities 
possible explanations are: 

Recent burstiness: The different level of bursti- 
ness is not surprising because these fields are sep- 
arated by a distance (850 pc) larger than most of 
the star forming complex es discovered in the SMC 
(see Livanou et al. l2007t ) , thus allowing the recent 
activity to fiuctuate independently; 

Recent systematic behavior: The stellar density 
in the SFH4 region is lower than in SFHl. Hence 
the systematic decrease of SF activity in SFH4 
may be connected with a minor amount of fuel 
available to support it up to today; 

Early quiescence and prompt rise: Are the qui- 
escent period and rapid SF increase 5-6 Gyr ago 
a global property of the Bar? Recent observations 
of RGB stars have revealed that older stellar com- 
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Fig. 20. — Bologna SFHs (black line) compared to Cole's (red line). The left panel and the right panel refer 
to the fields SFHl and SFH4, respectively. The inset panels zoom in on the last 1 Gyr. 



ponents of the SM C have a velocity dispers ion of 
about 27.5 km s"^ ([Harris fc Zaritskvlliooeh . high 
enough to distribute the stars over a large distance 
(of the order of few kpc) from their birth places 
within few Gyr. Hence, the low early activity and 
the prompt rise are not peculiarities of our fields 
but global features of the SMC Bar . Moreover, 
Subramanian &: SubramaniamI (|2012l ) find no ev- 
idence for a Bar in older stars, which is consis- 
tent with our low early activity. Further su pport 
is also provided bv lChiosi fc Vallenaril (|2007l ). who 
found similar star formation trends for three other 
Bar fields located around the SMC clusters K 29, 
NGC 290, and NGC 265. From the theoretical 
point of view it is not clear what mechanism is 
responsible for the rapid rise of stellar produc- 
tion. Was it externally triggered by the MW or 
the LMC, or self-initiated? The striking similar- 
ity between the SMC and LMC SFH is suggestive 
of the former. Pointin g in this direc t ion ar e the 
recent calculations by Diaz fc Bekkil (2011) who 
presented evidence that around 5.5 Gyr ago the 
LMC and SMC were within 160 kpc of the MW 
and 200 kpc of each other, therefore arguing for an 
independent origin of the Clouds. In this scenario 
the transition between quiescent and active phases 
could be naturally explained in terms of growing 



rate of mutual interactions started around 5 Gyr 
ago. 

Finally, iRafelski fc Zaritskvl (|2005l) noted that 
the spatial distribution of the younger and older 
clusters in the SMC is statistically different, lead- 
ing to the inference that a significant accretion or 
merger event may have taken place around 3—5 
Gyr ago. 

The study reported in this paper was the first 
step in a wider research activity aimed to charac- 
terize the SMC SFH through deep HST/ACS ob- 
servations. A forthcoming paper will be dedicated 
to the analysis of the Wing and Halo fields. This 
will allow a comparative analysis to look for global 
physical characteristics in the SMC star formation 
process. 
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Fig. 21. — Best-fit age-metallicity relation of the two Bar fields derived by Cole's method. The two fields show 
similar chemical evolution; within the errorbars they are the same. Taken together they paint a consistent 
picture of average metallicity increasing steadily over time. 
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